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Abstract 

In [19], an implementable algorithm was introduced to compute discrete solutions of 
sweeping processes (i.e. specific first order differential inclusions) . The convergence of this 
numerical scheme was proved thanks to compactness arguments. Here we establish that this 
algorithm is of order \ . The considered sweeping process involves a set- valued map given by 
a finite number of inequality constraints. The proof rests on a metric qualification condition 
between the sets associated to each constraint. 

Key-words: Differential inclusions - Subdifferential calculus - Numerical analysis - Prediction- 
correction algorithm. 
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1 Introduction 

In this paper, we are interested in the discretization of some first order differential inclusions 

§(t) + l7(f,q(t))9f(*,q(f)), (1) 
involving a multivalued operator U and a mapping f . 
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In the case where U is the sub differential of the inelicatrix of a closed convex set K (U(q(t)) = 
<9Ift;(q(i))), G. Lippold has shown in [11] that the implicit Euler method converges with an order 
h, where h is the time-step. In the framework of sweeping processes by moving convex sets K{t) 
(i.e. U(t, q(t)) = 9Ij^f t \(q(t))), J.J. Moreau has proved in [15] the existence of solutions of 
(PQ) with / = 0. More precisely, he introduced the Catching-up algorithm to construct discrete 
solutions converging with an order of precision ^. In [I], U takes the form of a maximal monotone 
operator perturbed by a Lipschitz coercive function. Similarly the authors regain that the order 
of convergence is ^. In [3], U is the proximal normal cone to a uniform prox-regular set Q(t) (i.e. 
U (t, q(i)) = N(Q(i), q(t)) (see Definition 12. 1|) ). Even if U is not maximal monotone, G. Colombo 
and V.V. Goncharov also recover the same convergence rate (with / = 0). These schemes 
adapted to differential inclusions have given rise to other works, see for instance [H [TH| [TUl 15] IS] . 

Here, we deal with a perturbed sweeping process: 

^(t)+N(g(t),q(t))9f(t,q(t)), (2) 

where Q(t) is defined by a finite number of inequality constraints and known to be uniform 
prox-regular. We recall that the projection onto such a set is well-defined in its neighbourhood 
but it cannot be exactly computed. That is why we consider a modified numerical scheme, that 
is based on a local approximation of Q(t). Near a point q G Q(t) the set Q(t) is replaced with a 
closed convex set Q(t, q). This substitution makes loose the prescribed regularity of the moving 
sets (actually these sets are generally supposed to vary in a absolutely continuous way in the 
sense of the Hausdorff distance). In order to go around this problem, we need to check metric 
qualification conditions between the sets associated to each constraint. In this context, we keep 
the same order of convergence: 

Theorem 1.1. There exists a constant C$ > such that for h small enough 

Wih - q|U°°([o,T]) < Coh 1 / 2 , 

where q and qh are the continuous and discrete solutions of (dj). 

We emphasize that this new approach allows us to avoid resorting to compactness arguments, 
used in |19| . So it permits to extend the convergence result of [19] in an infinite dimensional 
Hilbert space. 

The paper is structured as follows: In Section [21 we describe the mathematical framework by 
specifying notations and assumptions which will be used throughout the paper. Then in Sec- 
tion [3l after recalling the prediction-correction scheme proposed in [19], we prove in Theorem [33] 
that the discrete solution converges to the exact solution with order ^. This proof rests on a 
metric qualification condition which is checked in Section [H Finally, we illustrate this result 
with numerical simulations. 



2 Context 

In the sequel, the space M. d is equipped with its Hilbertian structure. We write B(x,r) for the 
open ball of center x G R rf and radius r > 0. 

We consider perturbed sweeping process by a set-valued map Q : [0, T] R d satisfying that 
for every t G [0, T], Q(t) is the intersection of complements of smooth convex sets. Let us first 
specify the set-valued map Q. For i e {1, • • • ,p}, let g t : [0, T] x R d -> R be a convex function 
with respect to the second variable. For every t G [0, T], we introduce the sets Qi(t) defined by 

Qi(t) := {q£R d , 9i(t,q) > o} , (3) 
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and the feasible set Q(t) (supposed to be nonempty) is 

Q(t) := flQiit). (4) 

i=l 

Let f : [0, T] x R d -> M d be a map, the associated (perturbed) sweeping process can be expressed 
as follows: 

^(t) + N(Q(t), q(t)) 9 f(t, q(t)) for a.e. t G [0, T] 
q(0) =q oG Q(0). 

We write N(Q(t),q(t)) for the proximal normal cone to Q(t) at q(i), defined below. 

Definition 2.1 Q2J). Lei 5 6e a closed subset ofR d . 
The proximal normal cone to S at x is defined by: 

N(S, x) := |« G M d , 3a > 0, s € P s (x + ot)J, 

Ps(y) : = G -5, <%(y) = \y - z\}, with d s {y) := inf |y - z| 
corresponds to the Euclidean projection onto S. 

We now come to the notion of uniformly prox-regular set. It was initially introduced by H. 
Federer (in [8]) in finite dimensional spaces under the name of "positively reached set". Then 
it was extended in infinite dimensional spaces and studied by F.H. Clarke, R.J. Stern and 
P.R. Wolenski in [2] and by R.A. Poliquin, R.T. Rockafellar and L. Thibault in [TB] . 

Definition 2.2. Let S be a closed subset ofM. d and rj > 0. The set S is said n -prox-regular if 
for allxeS and v G N(5, z) \ {0} 

B fx + ri^,7] \ n S = 0. 

Equivalently (see ]16p. S is r] -prox-regular if and only if the projection operator P5 is single- 
valued on 

{z e R d , d s (z) < r,}. 

This differential inclusion can be thought as follows: the point q(t), submitted to the pertur- 
bation {(t,q(t)), has to stay in the feasible set Q(t). To obtain well-posedness results for ([5]), 
we will make the following assumptions which ensure the uniform prox-regularity of Q(t) for all 
t G [0, T]. We suppose there exist c > and open sets Ui(t) D Qi(t) for all t in [0,T] such that 

d H (Qi(t),R d \Ui(t))>c, (AO) 

where dn denotes the Hausdorff distance. We set U(t) := Hf=i Ui(t). Moreover we assume there 
exist constants a, (3, M > such that for all t in [0, T], gi(t, •) belongs to C 2 (Ui(t)) and satisfies 

Vqet^t), a< \V qgi (t,q)\ < {3, (Al) 

VqGR d , \d tgi (t,q)\ <0, (A2) 

VqG^(t), \d t V q9i (t,q)\ <M, (A3) 

and 

VqeUi(t), \Dl gi (t,q)\ < M. (A4) 
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For all t G [0, T] and q G Q(i), we denote by q) the active set at q 

I(t,q) :={i€{l, ...,p}, 5i (t,q) = 0} 

and for every p > 0, we put: 

^(*,q) := {« e {i< •••,£>}, 3i(*,q) < p} • 
In addition we assume there exist 7 > 1 and p > such that for all t G [0, T], 



VqeQ(t), VAi >0, ^ Ai|V q ft(t,q)| < 7 

ie/ p (i,q) 

We will use the following weaker assumption: 



VqG Q(t), V\i > 0, ^ Ai|V q5 i(*,q)| <7 

ie/(t,q) 



X] AiV q 3i(£,q) 
*e/ P (t,q) 



A iV q 5i(t,q) 

ie/(t,q) 



(6) 
(7) 

(A5) 



(A5' 



In particular, this last assumption implies that for all i, the gradients of the active inequality 
constraints V q <7j(£,q) are positive-linearly independent at all q G Q(t), which is usually called 
the Mangasarian-Fromowitz constraint qualification (MFCQ). Conversely the MFCQ condition 
at a point q yields a local version of Inequality (|A5'|) . 

We recall some useful results established in [19] (Propositions 2.8, 2.9, 2.11 and Theorem 2.12 

in rrnn. 



Proposition 2.3. For all t G [0,T] and q G Q(t), 



N(Q(t),q) = ^N(Q i (t),q) = - £ R+V q5l (t,q). 

iel(t,q) 



Proposition 2.4. Under assumptions &A1\) . \A4\) and &A5]) . the set Q(t) is rj-prox-regular with 



n 







M7 



, for every t G [0, T] . 



Proposition 2.5. Under assumptions HA0\) . I(A1\) . HASty and HA5\) . the set-valued map Q is 
Lipschitz continuous with respect to the Hausdorff distance. More precisely there exists Kl > 
such that 

Vt,sG [0,T], d H (Q{t),Q(s)) <K L \t-s\. 
Theorem 2.6. Let T > and f : [0, Tjxl^ R d be a measurable map satisfying: 

3K f >0,Vq,qG (J Q{s) , Vt G [0, T] , |f(i, q) - f(t, q)| < K { \q - q| (8) 

se[o,T] 

3L f >0,VqG |J Q(s), Vt G [0,T], |f(i,q)| < L f (l + |q|). (9) 
s e[o,T] 

Then, under Assumptions iA0\) . \A1\) . HASfy . \A4\) and (Ad]) for all qo G Q(0), the following 
problem 

^(t)+N(Q(t),q(t)) 3 f(i,q(i)) /or a.e. t G [0,T] 
q(o) = qo, 

has one and only one absolutely continuous solution q satisfying q(t) G Q(t) for every t G [0,T]. 
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3 Time- stepping scheme 

Let us detail the numerical scheme proposed in [19] to approximate the solution of (jlOp on the 
time interval [0, T]. Let re G N*, h = T/n be the time step and t k = kh be the computational 
times. We denote by q k the approximation of q(t^) with qQ = qo. The next configuration is 
computed as follows: 

q2+i:=p^ liqB )(q* + «(«.q2)) ( n ) 

with 

p 

Q(t,q) :={qeR d , 9i (t,q) + (V q <7i(i,q),q- q) > V*} for q € I7(t) := f| Z7 4 (t). 

i=l 

We recall that all the gradients V q pi(t,q) are well-defined provided that q G U(t). Indeed it 
can be checked that this scheme is well-defined, more precisely for h < 

Q(t n k+1 ,q n k )cQ(t n k+1 )cU(t n k+2 ) (12) 

with c and Kl respectively given by Assumption ()A0p and Proposition 12.51 (see Proposition 3.1 
in |19j). Thus every computed configuration is feasible and the set Q(t, q) can be seen as an inner 
convex approximation of Q(t) with respect to q. The dependence on q of Q(t, q) do not allow us 
to have recourse to usual techniques (see for instance [151 Q2J E] ) based on the time-regularity 
of the sets Q(t). 

This scheme is a prediction-correction algorithm: predicted position vector q k + M(t k , q k ), that 
may not be admissible, is projected onto the approximate set of feasible configurations. 

Before stating the result of convergence, we introduce some last notations. We define the 
piecewise constant function f n as follows, 

P(i) = f(t n k ,q n k ) if t e k < n and P(T) = f(t^i,<&-x)- (13) 

We denote by q n the continuous, piecewise linear function satisfying for k G {0, . . . , re}, q n (t k ) = 
q k . To finish, we introduce the functions p and 6 defined by 

p n (t) = t n k and 9 n (t) = t n k+1 if t G [tl t n k+1 [, p n (T) = T and 6 n (T) = T. 

We recall some results about these approximate solutions (see Subsection 3.2 in [19] for details) : 

Theorem 3.1. Let us suppose that there exists Hf > for all q G (^J Q(s) and t,s G [0,T] 

se[o,T\ 

\i(t,q)-i(s,q)\<H i \t-s\ 1 / 2 . (14) 

Then with the assumptions of Theorem \2.(A q n tends to q in C°([0, T], W d ), where t h-> q(t) is 
the unique solution of U0\) . 

Remark 3.2. We can replace the definition H13\) of f n with 

f n (t) = T [ tk+1 f(s,q%)ds if te[t n k ,t n k+1 [, k<n and f n (T) = \ C ffoq^Jds. 

to Jt" to JT-h 

k 

Proposition 3.3. There exist C, D, K > such that 

dq n 



sup||q n ||Loo ( [ 0iT]) < C, sup 

n n 

and for n large enough, 



dt 



< K 

i°°([0,T]) 



d Q(tl +1 ,^) < Dh. (15) 
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We now come to the main result of the present paper which specifies the convergence order of 
the previous scheme. 

Theorem 3.4. Under there exists a constant Cq > such that for n large enough 

1/2 



|q n - q|U«([o,T|) - c °\n) ' 



where q is the solution of f5j). 



Proof. We check that the sequence (q n ) n is of Cauchy type. 
Let m > n be large enough. Since for k G {0, .., n — 1} 

and 0(*fc+i) c lfc) is a closed convex set, it comes : for all y G R d 

(<£ + - q2 +a ,y- q£ +1 > < K + hf(%,<$) - qj^l^^O/). (16) 

By Assumption © and Proposition 13. 3| we get for almost t G [0, T], 



do n 



< L f (l + C) + if. (17) 



Consequently by dividing (|16p by /i, we obtain for all y G M rf 



— (t) - P(*),y - q"(0»(t)) > < Cl^(^ ( t), q n (p n (t)) (2/) (18) 



with Ci := £/(H- C7) + 
Taking y = q m {9 m (t)), it follows 



fw-mmtFiTW)) < Ci^ w , q "(^( t) )(q m (^ m W))- 



First case : |q m (# m (t)) — q n (p"(i))| < r/8 (with r later introduced in Theorem 14. 9p . 
Let us choose w G P Q(e n (t) )(q m (6' m (t))). Hence w G Q(9 n (t)) and 

|« - q"(^(t))| < k - q m (^ m (t))| + ~ < d H (Q(nt)), Q(9 m (t))) + r -<^ + r -< r - 

by Proposition 12.51 for n large enough. Moreover from (|12|) it comes q n (p n (t)) G Q(p n (t)) C 
U(9 n (t)) for n > TKl/c and then Inequality f)15|) implies that 

r 

^Q(6»"(t),q"(p™(t))(q n (/ 9n (*))) - 4 

for n large enough. Then, by Theorem 14.91 and Proposition 13.61 we deduce 



d 



Q(6" itU n ipn(t)) H < K\q n (p n (t))-w 
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where k := QpM/(2a). Hence with Propositions 13.3 1 and [27 



d 



Q(en m n {p n m ^ m (0 m (t))) < K\0?{ P n (t)) - W \ 2 + W m (9 m (t)) - W 



< K\q n (p n (t)) - w\ 2 + d H (Q{9 n {t)),Q{9 m {t))). 
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Since 

|q>"(*)) - »l < |q"(c"(t)) - q"WI + |q"W - q m WI + |q°*W - q m (<H*))l + |q m (e m (t)) - M 
<^ + | q » (4) _ qm(t) l + ^ + ^ 

n mm 

<\q n (.t)-q m (t)\+2 + -^, 

n n 

we get 

d Q(^w, q "(p"W))( c i m ( 0m ( t ))) ^ 2 *li B (*) " <n*)l 2 + 2k ( 2KT+ n KLT y + ML. 



Finally, 



^(t) - P(i), q m (9 m (t)) - q n (6 n (t))) < 2C lK \ q n (t) - q m (t)\ 2 + % 
at n 



with C 2 := Ci(K L T + 2kT 2 (2K + K L ) 2 ). 

Second case : |q m (6» m (t)) - q n (p n (t))\ > r/8. 
Then by (TBI). 



^(i) - r(i),q ro (r(0) - q"(0" < Ci|q m (0 m (t)) - q n (^(t))| 

< d|q m (^(t)) - q"(p n (*))l + d — 

n 

< ^C 1 |q m (0 m (t)) - q n (p n (t))\ 2 + C X — 
r n 

16 „ , _ n 2 64 „ /i<CT\ 2 „ 



<-C 1 \ q m (t)- q n (t)\ z + -C l [ +d 

r r \ n J n 

End of the proof : 

By setting C 3 := max{2CiK, 16Ci/r} and C 4 := max{C 2 , 64Cl ^ T ) 2 + Ci-fiTT}, we get 

^(t) - P(i),q m (# m (t)) - q"(0"(t))\ < C 3 \ q m (t) - q n (t)\ 2 + ^. 
Therefore 

^(i) - q m (t) - q»(t)\ < C 3 |q m (t) - q n (t)| 2 + ^, 

with C5 := C4 + 2C\KT. By summing the previous inequality and the other one obtained by 
changing the role of n and m, it yields 

' (*) - ~ir(t), q m (t) - q n (t) ) < 2C 3 \ q m (t) - q n (t)\ 2 + + |P(i) - f m (t)\\ q m (t) - q n (t)\. 



dt dt I n 

Furthermore by ([8]) and ()14[) 

|f"(t) -f m (t)| < \{(p m (t), q m (p m (t))) -f(p m (t),q"(p"(t)))| + |f(p m (t),q"(p"(t))) -f(p"(t),q n (p n (t)))l 

< ^ f |q m (p m (t)) - q n (p n (t))\ + H { \p n (t) - p m (t)\ l 2 

< K { \ q m (t) - q n (t)\ + 2K^- + H f \p n (t) - p m (t)\K 
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and so 

|P(t) - f m (t)||q m (t) - q n (t)| < K f \q m (t) - q n (t)\ 2 + + H t (^\ \q m (t) - q n (t)\ 



< (~ + ^ f )|q m (t) - q"(t)| 2 + ~ (2K{KT + H.Vf) 2 ~. 



1 (t)-^(t),q™(t)-q«(t)\<(2C73 + TT+^f)|q m (t)-q n (t)| 2 + ' 



Finally, 

with C 6 := 2C 5 + \ (2K{KT + HiVfY. 
By applying Gronwall's Lemma, we have 

1 /2 

l|q m - q n |U-([o ) T]) < (^) exp((2C 3 + ~ + K { )T). 
Then, we conclude the proof by taking the limit for m — > oo. □ 

Remark 3.5. This proof allows us to get around the compactness arguments employed in 119)/ 
to obtain the convergence ofqt- Consequently, this result can be extended to the Hilbertian case. 
Then it can be checked that the limit satisfies the differential inclusion ([3JJ by following the same 
reasoning as in [19]. 

It remains to prove Proposition 13.61 and Theorem 14, 91 We now check the first result whereas the 
second one will be established in the next section. 

Proposition 3.6. Fori G {l,...,p} and q G U{t) we set 

Qi(t,q) :={qeR d , 9i (t, q) + (V q 9i (t, q), q - q) > 0}. 

Then, for all t G [0,T], all q G Q{t) and all q G U(t), we have for all i G {1, ■■■,p} 

^ ( ^(q)<^|q-q| 2 - (is) 

Proof. Let consider i G {1, •••,£>}, q G U(t) and q G Q(t) C Qi(t). We assume that q G" Qi(t, q) 
(otherwise (fTUj) obviously holds). 
For I > 0, we define 

Z(£) :=q + eV gi (t,q). 
The point z(£) belongs to Qi(t,q) if and only if 

9i {t, q) + {V 9i (t, q), q - q) + ^|Vffi(t, q)| 2 > 0, 

which is equivalent to 

9i (t,q) + (V 9i (t,q),q - q) 



Thus, 



dQ dm) (q) < |q- z(£ )\ <£o|Vft(t,q)| 

3i(*,q) + (V^(t,q),q - q) 



< 



|V gi (t,q)| 



1 

—-7 / sDq5i(t,q + s(q-q))(q-q,q-q)<is, 
because 9 i(t,q) > 0. We conclude to ([THD by Assumptions (|Aip and (|A4|) . □ 
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4 Metric qualification condition 

This section is devoted to the proof of Theorem 14,91 which corresponds to a metric qualification 
condition for the sets Qi. Aiming that, we recall some notions of subdifferential calculus. 

Definition 4.1 (proximal subdifferential). Let f : M. d — > R be a lower semicontinuous function 
which is finite at x G M. d . The proximal subdifferential of f at x is defined by 

8 p f(x) := jx* G R d , 3a, p > 0, V\h\ < (3, f(x + h) - f{x) > (x*, h) - a|/i| 2 } . 

Definition 4.2 (limiting subdifferential). Let f : W 1 — > R be a lower semicontinuous function 
which is finite at x G M d . The limiting (or Mordukhovich) subdifferential of f at x is defined by 

d L f(x) := I x* G R d , x* = lim x* k with x* k £ d P f(xk), Xk -)■ x and f(x k ) -)■ f(x) \ . 

y k— >oo ' J 

Definition 4.3 (Clarke subdifferential). Let f : M. d — > R be a Lipschitz continuous function. 
The Clarke subdifferential d f(x) of f at x can be defined (see J^) as the closed convex hull of 
the limiting subdifferential : 

d f(x) := conv d L f(x). 

This notion has been extended for less regular functions, we refer the reader to [T7j for details. 

The following property is a special case of the exact sum rule for the Clarke subdifferential (see 
Theorem 2 of [17J): 

Lemma 4.4 (Optimality property). Let f : M. d — > MU{+oo} be a lower semicontinuous function 
and (f) : W 1 — > M a convex Lipschitz function. If x £ M rf is a finite local minimum of f + (/> then 

G d c f(x) + d c 0(x). 
Let us recall the variational principle of Ekeland (see [7]). 

Proposition 4.5 (Ekeland variational principle). Let f : M, d — > M U {+00} be a lower semi- 
continuous function which is bounded from below. Let e > and x G M. d such that 

inf / < f{x) < inf / + e. 

Then for all A > ; there exists w G M. d satisfying 

• f(w) < f{x) 

• \x — w\ < A 

• for all z / w, f(z) > f(w) — j\z — w\. 

The following result comes from Theorem 2.1 in [SJ. For an easy reference, we detail the proof. 

Lemma 4.6. Let f : M. d — > M + U {+00} be a lower semi- continuous function and xq with 
f(xo) = 0. Assume there exist 7, 5 > such that for all 

x * G |J d C f(x) 

f(x)>0 

we have \x*\ > 7. Then for all x G B(xq, 5), d{f =0 }(x) < 7 _1 /(x). 
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Proof. Let x £ B(xo,5). If /(x) > 7<5, then 

d{ f=0 y(x) < \x - x | < 6 < ^ l f(x). 

Now, we assume that < f(x) < j5 and we set e := /(x). Applying the variational principle 
of Ekeland (see Proposition I4.5P to / with e and any A £]7 _1 e,<5[. There exists w = w(X) £ M d 
such that f(w) < f(x), \x — w\ < A and 

Vz/u>, f(z)>f(w)-j\z-w\. 

Consequently, w minimizes / + eA _1 | • —w\ and by Lemma 14.41 it comes 

£ d c f{w) + d c ?p(w) 

where ip(-) = eA _1 | • —w\. 

So there exists x* £ d f(w) with |x*| < eA _1 < 7. If f(w) > 0, that is in contradiction with the 
assumptions as \w — xq\ < \w — x\ + \x — xo\ < 26 and so we deduce that necessarily f(w) = 0. 
Then we conclude to the desired result, since 

dff=o}(x) < \x — w\ < A 

holds for every A £]7 e, 6[. □ 

From now on, we come back to the framework of the previous sections and prove the metric 
qualification condition of sets Qi . 

v 

In the sequel, we introduce convex sets Cj for i £ {1, ...,p} and their intersection C = f\Ci. 

i=i 

We consider the following set-valued map F 

\x ^ := (Ci-x) x ••• x (C p -x). { ' 

Let us note that £ F(x) if and only if x £ C. 

Proposition 4.7. Consider the function f defined by f(x) := rfjr(a;)(0) where F is given by I12(J\) . 
The map f is Lipschitz continuous and for all x ^ C, 



d P f(x)cd C f(x)={ £ |- 

2 5 cc(^. C j 



where y = Ppf x \(Q). In other words, for all i £ {1, ...,£>}, yi + x £ Pc^x), hence —yi £ 
N(Ci,x + w)- 

Proof. For all x £ M d , 

/(x) = d F{x) (0) = d n (4>(x)) 
where II := <8>f =1 Cj and </>(x) := (x, . . . , x) £ For x ^ C, 

d c f(x) = d c (d n o 0)(x) = . . . , 1) • 5 c d n (0(x)) 

thanks to Corollary 1 in |17j . By convexity of the sets Ci, du is a convex function and so 

d c f(x) = t (l,...,l)-d p du(Hx)), 
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see Remark 14.81 First we claim that 

d p d u (Hx)) C [g^^C^ri 5(0,1), (21) 

with £i(x) := -^^)d p d Ci (x) \i x ^ d and £j(a:) := {0} else. 

Indeed, let x* belong to d p dn{4>(x)) . By definition, for some a > and for all small enough 
h G 

dn(<^(x) + /i) - d n (0(x)) > (x*, M - a|/i| 2 . 

Let us fix an index i G {1, . . . ,p} and set h = (0, . . . , 0, hi, 0, . . . , 0). Since dn(</>(")) 2 = Si > 
it follows that for all small enough hi G M rf 



^d u ((f>(x)) 2 + d Ci {x + /ii) 2 - d Ct {x) 2 - ^d n (4>{x)) 2 > {x*, hi) - a\hi\ 2 . 

By a first order expansion, we get 

d Ci (x + hj) 2 - dciixf , * , , _ , 2 
2d n (^)) ^^'^ Q£|/l » l ' 

with another numerical constant a' . Then, we obtain with another constant a" and for all small 
enough hi G M d 

(dCiOc + /ii) - dc^x)) > (x*, hi) - a"\hi\ . 



du{<t>(x)) 

If x G Ci then (x) = and so we deduce that x* = 0. If x ^ Cj then by definition of the 
proximal normal cone, 

xt£d P d Cl (x)cS(0,l), 



X Zi 



dCi (x) 
see Remark l4~8l 

So \x*\ = dc i {x)dn{4 , {x))~ 1 and so \x*\ = 1, which concludes the proof of ([2T]) . 
Let us now finish the proof of the proposition. Thus 

We set 2 = (zi,...,,Zp) G with for all i, Zj = Pc.(i) or equivalently z = Pn(0(x)). By 
Theorem 1.105 in |14| . 

d p d Ct (x) cd p dcM)n 5(o,i) f ' l ' 

Consequently, we have 

We finish the proof by choosing y := <p(x) — z G M. dp . □ 

Remark 4.8. Let S C M. d be a closed convex set and x ^ S, then d p ds{x) = d c d$(x) C 
5(0, 1). Indeed with w := Ps( x ) and vectors h = e(w — x) for small enough e, we remark that 
ds(x + e(w — x)) = dg(x) — e\w — x\ . Hence, by Definition \4-l\ we obtain for every x* G d p ds{x) 

-\h\ > (x*,h) -a\h\ 2 . 

By dividing by \h\ and letting e go to 0, we deduce that \x*\ > 1. We also conclude to \x*\ = 1 
since ds is 1-Lipschitz. 
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Theorem 4.9. There exist r and 6 such that for all t G [0, T], for all q 6 U{t) satisfying 
d Q(t,q)^ - r / 4 and al1 q G - B (^' r / 4 )' 



Indeed we can choose O = and r = min( 



i=l 



4p g 



13/3 ' 2M7 / ' 



Proof. Consider r = min(^, jm^)- Let us fix i G [0,T] and q G satisfying dg^-^q) < r/4. 

Consequently there exists qi G i?(q, r/4) such that qi G q). We define a Lipschitz map 

/ := dp(.)(0) where F is given by ([2TI|) with C, = q). First we check the assumptions of 

Lemma HTBl for the function / with xo = qi- Indeed /(qi) = because qi G Q(t,q). 

Let us consider q G B(qi,r) n Q(^q) C ; so q G B(q,2r). By Proposition 14. 7| d c f(q) = {q*} 

where 



q := 



«, q£<9i(*,q) 



with p* = p/|p| and p = P_F( q )(0). Moreover for i satisfying q ^ Qi(t,q), — p* G N(Cj,q + p« 
Let us define 

•J(*,q) :={j, 9j(t,q) + (V^-(t,q),q- q) < 0} = {j, q^Qj(t,q)}. 
It is well-known that there also exist nonnegative reals (Aj)j e j(£ jq ) satisfying 



q* = Wg t (t,q). 

i£j(t,q) 



Hence by Assumption (|A4 



> 



AiVftfoq) 

i€J(t,q) 

^ AiVft(t,qi] 

i£J(t,q) 



Mr x - 



ieJ(t,q) 

Since q + p» G Pc 4 (q), Proposition 13,61 yields |pj| = dc;(q) < f^(2r) 2 < r. Moreover for all 
i G J(i, q), q + Pi G dQi(t, q) so we have by Assumption (|Aip with the convexity of ^ 

Si(i,q) < -(Vc/i(i,q),q + Pi - q) < /3(|q - q| + |Pi|) 
< 3/3r. 
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Hence by Assumption (jAljl . 



Si(*,qi) < 3/3r + /3r/4 = — /3r 



Due to the choice of r, we deduce that 13/3r < \p and thus J(i, q) C I p (t, qi). From Assumptions 
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TJ), (JA4I) and (|X5j) . we deduce that 



|q*|>7 _1 Yl ^IVft^qi)!-^ ^ Ai 

i6J(*,q) «eJ(t,q) 



i£j(t,q) 



^-^> E * 



i€J(t,q) 

^ E Ai|V«(t,q)| 

,P ieJ(t,q) 

a \ ^ . .. 

^ E ip*i 

»6J(t,q) 

a . a 



27/3" 1 27/3 

We can also apply Lemma 14.61 and we obtain that for all q € £?(qi,r/2) D i?(q, r/4) 

(\ 1/2 
ieJ(t,q) I i=l 



□ 



5 Numerical simulations 

The aim of this section is to illustrate the convergence order with an example (due to the 
modelling of crowd motion in emergency evacuation). We refer the reader to [18} H2] for a 
complete and detailed description of this model. 

We quickly recall the model. It handles contacts, in order to deal with local interactions between 
people and to describe the whole dynamics of the pedestrian traffic. This microscopic model 
for crowd motion (where people are identified to rigid disks) rests on two principles. On the 
one hand, each individual has a spontaneous velocity that he would like to have in the absence 
of other people. On the other hand, the actual velocity must take into account congestion. 
Those two principles lead to define the actual velocity field as the Euclidean projection of 
the spontaneous velocity over the set of admissible velocities (regarding the non-overlapping 
constraints between disks). 

More precisely, we consider N persons identified to rigid disks. For convenience, the disks are 
supposed here to have the same radius r. The center of the i-th disk is denoted by q« G M 2 . 
Since overlapping is forbidden, the vector of positions q = (qi, ..,q7v) G M. 2N has to belong to 
the "set of feasible configurations" , defined by 

Q := {q G R 2N , Ai(q) > Vi ^ 3} , (22) 

where .Djj(q) = |q^ — qj\ — 2r is the signed distance between disks i and j. If the global 
spontaneous velocity of the crowd is denoted by U(q) = (U\(qi), .., Un(<In)) G M. 2N , the previous 
crowd motion model can be described by the following differential inclusion: 

§+N(Q,q)9U(q). (23) 
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This evolution problem fits into the theoretical framework developed in this paper. 

For the numerical simulations, we treat an emergency evacuation of a square room (10x10) 
initially containing N = 150 persons (identified to rigid disks of radius r = 0.2). Since the exact 
solution is unknown, we approach the error as follows 

Hq-qfclU- - max n Jqh min (ti)-qA(ti)l := e h (24) 
1=1,... ,10 

where t { := »T/10. We choose h min := 0.01 and h € {0.02, 0.025, 0.04, 0.05, 0.0625, 0.08, 0.1, 0.2, 0.5}. 
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Figure 1: Log-log plot of the error versus h. 

We observe in Figure Q] that the empirical order of convergence is \ ■ Note that for a time 
step h close to h m i n , the approximation (j24"|) is not valid. 
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